Non-Newtonian granular hydrodynamics. What do the inelastic simple shear flow and 

the elastic Fourier flow have in common? 
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We describe a special class of steady Couette flows in dilute granular gases admitting a non- 
Newtonian hydrodynamic description for strong dissipation. The class occurs when viscous heating 
exactly balances inelastic cooling, resulting in a uniform heat flux. It includes the Fourier flow 
of ordinary gases and the simple or uniform shear flow (USF) of granular gases as special cases. 
The rheological functions have the same values as in the USF and generalized thermal conductivity 
coefficients can be identified. These points are confirmed by molecular dynamics simulations, Monte 
Carlo simulations of the Boltzmann equation, and analytical results from Grad's 13-moment method. 

PACS numbers: 45.70.Mg, 47.50.-d, 51.10.+y, 05.20.Dd 



The study of granular matter is interesting from a tech- 
nological point view because its understanding has many 
applications in technology and in sciences other than 
physics, such as biology [HQ. Furthermore, it is also im- 
portant from a more fundamental point of view Q . For 
example, the Boltzmann equation (BE) for low-density 
granular gases describes a generalization of the contrac- 
tion from a microscopic to mcsoscopic scale and thus a 
generalization of fundamental concepts in the fields of 
statistical and fluid mechanics. The BE for granular 
gases (usually modeled as smooth inelastic hard spheres) 
has been widely employed to analyze several granular 
flow problems and a large number of research works have 
been recently published in this field jj H| . A standard 
approach used for solving the BE for ordinary gases con- 
sists in obtaining a perturbativc solution, which results in 
Navier-Stokes (NS) or Burnett type hydrodynamic equa- 
tions Q. However, the kinetic energy loss in the colli- 
sions renders thegranular steady flows inherently non- 
Newtonian 0, i, Hf . 

We report in this Letter on strong evidence of non- 
Newtonian hydrodynamic steady states in the planar 
Couette flow geometry for a wide range of inelasticities. 
Furthermore, our description is inclusive in the sense that 
it comprises a class (manifold) of steady flows, whose el- 
ements correspond both to granular and ordinary gases. 
This novel class occurs when the heat flux q is con- 
stant across the system, due to an exact local balance 
of inelastic cooling and viscous heating, even though the 
temperature and the shear rate are in general not uni- 
form. As a consequence, this class of nonlinear flows 
has several surprising hydrodynamic properties: (i) the 
nonlinear temperature profiles T(y) are indistinguishable 
from those of the Fourier steady state of an ordinary gas 
with the same temperature difference (Fourier flows are 
the "elastic-limit" elements of the manifold); (ii) when 
the spatial coordinate y normal to the moving plates is 
eliminated between temperature and flow velocity the re- 
sulting profiles T(u x ) are linear; (hi) the non-Newtonian 
rheological properties (shear stress and anisotropic tem- 



peratures) are uniform and have the same values as those 
obtained in the well known simple or uniform shear flow 
(USF) of granular gases 0, which is in fact the partic- 
ular case q = of the manifold; and (iv) the heat flux 
vector is proportional to the thermal gradient (general- 
ized Fourier's law) with an effective thermal conductivity 
tensor. Because of property (ii), we will refer henceforth 
to these flows as "linear T(u x ) flows," or simply, "LTu" 
flows. Property (iv) can be interpreted as a method for 
measuring the intrinsic thermal conductivity coefficients 
of the USF state directly from LTu steady states. All 
these results are supported by three independent and 
complementary routes: an approximate analytical solu- 
tion from Grad's 13-moment (G13) method to the BE, 
direct simulation Monte Carlo (DSMC) numerical solu- 
tions of the BE, and molecular dynamics (MD) simula- 
tions. The existence of the special LTu class at NS order 
was theoretically proven in a recent work Q, but the 
applicability of the NS description is restricted to the 
quasi-clastic limit and so the general proof of LTu states 
requires a non-Newtonian description, as carried out in 
this Letter. 

Conservation of momentum in the steady state Couette 
flow implies Pi y = const, where Pij is the stress tensor. 
As for the energy balance equation, it reads 



-d y q y = (3/2)nTC + P xy d y u x , 



(1) 



where n is the number density and £ is the inelastic cool- 
ing rate. Equation (JTJ is valid for all steady states in the 
system, whether hydrodynamics applies or not. As can 
be seen, the signature of the heat flux gradient is deter- 
mined by the balance between two terms: the first one 
on the right-hand side comes from inelastic cooling and 
is inherently positive; the second term is due to viscous 
heating and is inherently negative. Thus, the condition 
for homogeneous heat flux is that these two terms exactly 
balance in the whole bulk region. Moreover, the stream- 
wise heat flux component q x (absent at NS order) turns 
out to be "synchronized" to the crosswise component q y 
in the sense that it becomes homogeneous when so does 
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(Equilibrium) 

FIG. 1: (color online) Each point of this diagram represents 
a Couette flow steady state. The surface defines the LTu 
class, which contains the lines representing the Fourier flow 
for ordinary gases and the USF for granular gases. 



Application of the classical G13 method [8| to the BE 
for inelastic collisions yields £ = j^v{1 — a 2 ) (where 
a is the coefficient of normal restitution and v — 
^■ypKna 2 \JT /m is an effective collision frequency, a and 
m being the diameter and mass of a sphere, respectively) 
and a closed set of coupled equations for the hydrody- 
namic fields, the stress tensor, and the heat flux. This set 
allows for an LTu solution characterized by nT = const, 



v 1 d y T = A = const, 



-n(a)d y u x = const, 



P 

nT 



-X(a)d y T = const, q x 



a(a) = const, (2) 



i(a) = const, (3) 



(a)d y T = const. (4) 



Since the trace of the stress tensor is 3nT, one has 
X + y + Z = 3. Equation © implies that dT/du x = 
A/a{a), so that T is indeed a linear function of u x . Equa- 
tions ©"(111) define a class of solutions because the con- 
stant A is arbitrary. On the other hand, the dimension- 
less shear rate a(a) (which is the Knudsen number as- 
sociated with the shearing), the effective shear viscosity 
7y*(a) = 77(a)/?7Ns(l), the temperature ratios Oi(a), and 
the effective thermal conductivities A* (a) = A(q;)/Ans(1) 
and 4>*(a) = 4>{a)/ Ans(1) ar e independent of A, their 
values depending on inelasticity only. Here ?7ns(1) = 
nT/v and Ans(1) = 15?7Ns(l)/4m are the NS transport 
coefficients of the elastic gas [8[. The quantities a(a) and 
T]*(a) are related each other by the exact balance equa- 
tion ([T]) and the LTu condition d y q y = 0, 



a 2 (a) = 3(*(a)/2r]*(a) 



(5) 



where (*(a) = C,/v. Equation ([5]) shows that in the elas- 
tic limit a — > 1 the LTu shear rate vanishes (a —> 0) and 
thus the conventional Fourier flow (4 ^ 0) for an ordi- 
nary gas is included in the LTu class as a special case. 
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FIG. 2: (a) Profiles of the heat flux component q y for a = 0.7, 
AT — 5, and a shear rate 8% smaller than the threshold 
value (a), equal to the threshold value (□), and 4% larger 
than the threshold value (■). (b) Temperature profiles in the 
LTu flow for a common fluid temperature difference T(h/2) — 
T(-h/2) = 4 and ct = 0.5 (□), a = 0.7 (A), and a = I (a). 
In the latter case (ordinary gas), the simulated state is the 
conventional Fourier flow (without shearing) . The inset shows 
T 3/2 (y). (c) T(ux) profiles for the LTu class with T(h/2) - 
T(-h/2) = 4 and a = 0.5 (□) and a = 0.7 (A). The shown 
data have been obtained by DSMC simulations. 



Conversely, the USF is recovered as another special case 
in the limit A — > with a < 1. This is sketched in Fig. 
[1] The explicit expressions predicted by the G13 method 
for the transport coefficients are Q 



0i 



Pi + 3C 



A* = /3 2 (706» y - 20 + 63C7/3 2 )/(50/3f - 63a 2 ), 



(j>* = 7a(210 y - 6 + 10/3 2 r/*)/(50/3 2 - 63a 2 ), 
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where ft (a) = (l + a)(2 + a)/6 and ft(a) = (l + a)(49- 
33a)/32. It must be noticed that rf(a) ^ V^ksi ) — 
(ft + ±C)- 1 and A*(a) + X* m (a) = (ft - 3D" 1 . 

To validate the G13 theoretical predictions, we have 
performed DSMC simulations of the BE and MD sim- 
ulations (global solid fraction = 7 x 10~ 3 ) for a granu- 
lar gas of hard spheres enclosed between two plates lo- 
cated at y = ±h/2 and moving with velocities U± (see 
Ref. Q for technical details). Diffuse boundary condi- 
tions characterized by wall temperatures T± (T_ < T + ) 
are applied. In what follows, quantities arc nondimcn- 
sionalized by the choice of units m = 1, T(—h/2) = 1, 
n(—h/2) = 1, and v(—h/2) = 1. In these units, the quan- 
tity A represents the maximum value across the system 
of the Knudsen number associated with the thermal gra- 
dient 0] . The separation between the plates has typically 
been set h ~ 5-20 and we have considered a wall tem- 
perature difference in the range AT = T + /T_ — 1 = 
20. We have looked for Couette flows belonging to the 
LTu class by fixing AT and varying the applied shear 
7 = (U + — U-)/h. Once the steady state is reached, we 
monitor the parametric plot of temperature versus flow 
field, T(u x ). We have observed in all the cases a definite 
sign of the curvature parameter d 2 T/du 2 in the bulk, 
with no inflection point. Interestingly, as the shearing in- 
creases and a certain threshold value 7 = 7th is crossed, 
the sign of d 2 T/du 2 undergoes a change from negative 
to positive. At the threshold shear 7th, d 2 T/du 2 x = 
and this signals the onset of the LTu flow, as explained 
above. This transition is accompanied by a change in 
the slope of q y , so one also has q y = const at 7 = 7th- 
This is illustrated by DSMC data in Fig. Oa), while 
Figs. [2Kb) and [U(c) show some representative LTu tem- 
perature profiles. Figure [2jb) is especially noteworthy 
since it clearly shows that all the LTu T(y)-profiles shar- 
ing the same temperature values near the walls collapse 
into a common curve independent of the inelasticity of 
the particles. Therefore, the temperature profile reached 
by the granular gas in the LTu flow is indistinguishable 
from that of an ordinary gas in the conventional Fourier 
flow. This surprising result is a consequence of the ap- 
plicability of hydrodynamics to granular gases, even with 
strong inelasticity. According to the first equation in ([2]) , 
d y T 3 / 2 = AvT 1 / 2 = const, so T 3 / 2 (y) is a linear func- 
tion that is completely fixed by the values near the walls, 
regardless of the value of a [see inset in[2Kc)]. On the 
other hand, since dT/du x — A/a(a), the slope in the 
T(it x )-profilcs is a-dependent, as shown in Fig.[2Kc). 

For each value of a we have computed the reduced 
shear rate a defined by the second equation of @, as 
well as the generalized transport coefficients defined by 
Eqs. ([3]) and (|4]). While the threshold value jth depends 
on AT, we have observed that, as predicted by theory, 
a(a), rf(a), 9i(a), A* (a), and <fi*{a) are insensitive to 
the choice of AT. In particular, as shown in Fig. [31 the 
reduced shear rate a(a) and the rheological quantities 
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FIG. 3: (color online) Plot of a(a), r/*(a), 9 x (a), and y (a) 
in the LTu flow as obtained from DSMC simulations (h = f 5) 
with AT = 2 (x) and AT = fO (+), and from MD simu- 
lations (h = 7) with AT = 2 (A) and AT = 5 (□). Also 
DSMC data (O) of the USF (AT = 0) QU are included. The 
lines represent the analytical results obtained from the G13 
approximation. Note that both in theory and simulation the 
values of 8 y axe hardly distinguishable from those of rf . 
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FIG. 4: (a) Plot of q x vs A for a — 0.7 (A) and of -q y vs A for 
a = 0.9 (a), as computed from DSMC simulations. The lines 
are linear fits, whose ordinate values at A = are 1.7 x 10 -3 
and 5.4 x 10 -3 , respectively, (b) Heat flux coefficients A* (a) 
(A, ■) and 4>*(a) (A, □) as obtained from DSMC simulations 
(triangles) and MD simulations (squares). The lines represent 
the analytical results obtained from the G13 approximation. 



if (a) and Oi(a) are the same in the LTu-class of Cou- 
ette flows (regardless of the value of AT) as in the USF 
(AT = 0) [IfJ, even though the boundary conditions 
are quite different: boundary-driven in the case of the 
Couette flo w lffl| and Lees-Edwards periodic boundary 
conditions [I2J in the USF case. Figure [3] also shows the 
close agreement between DSMC and MD results as well 
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as the reliability of the theoretical predictions from the 
G13 approximation. 

Now we turn to the heat flux coefficients. One of the 
most striking theoretical predictions is the linear rela- 
tionship between both components of the heat flux and 
the thermal gradient (generalized Fourier's law), as de- 
scribed by Eq. ((4]). This means that qi cx A and this is 
illustrated in Fig.HJa) for a = 0.7 and a = 0.9. The (re- 
duced) heat flux transport coefficients are evaluated as 
A* (a) = -{2m/5nT)q y /A and <jf{a) = (2m/5nT)q x /A, 
and are plotted in Fig. @Ib). It can be observed that 
the streamwise component q x becomes larger in magni- 
tude than the crosswise component q y for a < 0.9, what 
represents a strong non-Newtonian effect. Interestingly, 
this effect, as well as the general dependence of the trans- 
port coefficients are very well captured by our simple G13 
approximation. As happens with the rheological proper- 
ties, Fig. |4jb) shows a good agreement between DSMC 
and MD data for the generalized thermal conductivities. 

To sum up, we have described a special class of steady 
Couette flows (LTu class) in a low-density gas of inelastic 
hard spheres. This state encompasses the Fourier flow of 
clastic particles {d y T ^ 0, a = 1) and the USF of in- 
clastic particles (d y T — 0, a < 1) as special cases. In 
this sense, LTu can be seen as a "natural" extension (i) 
of the conventional Fourier flow in ordinary gases to the 
realm of granular gases and (ii) of the granular USF to 
states with nonzero heat flux. Therefore, the LTu uncov- 
ers a wide spectrum of inelasticities and wall tempera- 
tures within a unified framework, for both granular and 
elastic gases (see Fig.[T]). Three complementary and inde- 
pendent approaches have been followed: an approximate 
solution based on the Gf 3 method, DSMC simulations of 
the BE, and MD simulations of a dilute system. Here, in 
contrast to what happens in the Couette flow for ordinary 
gases [l^ . the G13 theoretical results compare surpris- 
ingly well with computer simulations, even for strong val- 
ues of dissipation [see Figs.[3]and[2Jb)]. This paradoxical 
result seems to be closely tied to the condition of uniform 
heat flux so that, as the balance between inelastic cooling 
and viscous heating breaks down, the Gf 3 approximation 
is not expected to give good quantitative results, even 
though a hydrodynamic description still applies. The 
solution found here clearly shows strong non-Newtonian 
effects since the shear viscosity and thermal conductivity 
coefficients qualitatively differ from their corresponding 
NS values and there exist anisotropic normal stresses and 
a nonzero streamwise heat flux. Despite the fact that 
q = in the USF, this state possesses intrinsic thermal 
conductivities in the same way as an ordinary fluid at 
equilibrium has intrinsic transport coefficients. Our re- 
sults on the LTu class assigns a meaning to the heat flux 
transport coefficients of the USF, which to our knowledge 
have been computed here for the first time. 

Since MD simulations just solve Newton's equation of 



motion for a many-particle system, the good agreement 
found here between the latter method and the DSMC 
method shows that the existence of the class of LTu flows 
reported in this Letter is not an artifact of the BE, which 
is based on the absence of spatial and velocity correla- 
tions (molecular chaos assumption). This fact can stim- 
ulate experiments with Couette geometry [lij ]. where it 
would be possible to test whether the linear relationship 
between temperature and flow velocity is attainable or 
not. Finally, it must be stressed that the analysis carried 
out here provides a nontrivial example of the existence 
of a hydrodynamic description for a strongly inhomoge- 
neous state beyond the NS regime for a dilute granular 
gas of hard spheres. 
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